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ABSTRACT 

Recent numerical relativistic simulations of black hole coalescence suggest 
that in certain alignments the emission of gravitational radiation can produce 
a kick of several thousand kilometers per second. This exceeds galactic escape 
speeds, hence unless there a mechanism to prevent this, one would expect many 
galaxies that had merged to be without a central black hole. Here we show that 
in most galactic mergers, torques from accreting gas suffice to align the orbit and 
spins of both black holes with the large-scale gas flow. Such a configuration has 
a maximum kick speed < 200 km s _1 , safely below galactic escape speeds. We 
predict, however, that in mergers of galaxies without much gas, the remnant will 
be kicked out several percent of the time. We also discuss other predictions of 
our scenario, including implications for jet alignment angles and X-type radio 
sources. 

Subject headings: black hole physics - galaxies: nuclei - gravitational waves - 
relativity 



Introduction 



When two black holes spiral together and coalesce, they emit gravitational radiation 
which in general possesses net linear momentum. This accelerates (i.e., "kicks") the coa- 
lescence remnant relative to the initial binary center of mass. Analytical calculations have 
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)), but because the majority of the kick is produced be- 



tween plunge and coalescence, fully general relativistic numerical simulations are necessary 
to determine the full recoil speed. 
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Fortunately, the last two years have seen rapid developments in numerical relativity. 
Kick speeds have been reported for non-spinning black holes with different ma ss ratios 
( [Herrmann. Shoemaker, fe Lagunal 120061 ; iBaker et al.l 120061 ; IGonzalez et al.l 120061) and for 
binaries with spin axes parallel or an tiparallel to the orbital axes (IHerrmann et al.l 120071 ; 



Koppitz et al. 



2007 



orientations ( Gonzalez et al. 



2007 



Baker et al.l 120071) . as well as in i tial ex plorations of more general spin 



Campanelli et al.l l2007al ibl) . For mergers with low spin 



or spins both aligned with the orbital angular momentum, these results indicate maximum 
kick speeds < 200 km s -1 . Remarkably, however, it has recently been shown that when 
the spin axes are oppositely directed and in the orbital plane, and the spin magnitudes are 
high (dimensionless angular moment um a = cJ/GM 2 ~ 1), the net kick speed can per- 
haps be as large as ~ 4 000 km s _1 ( IGonzalez et al.l 120071 : ISchnittman fe Buonannol 120071 : 
Campanelli et al.ll2007bh . 



The difficult y this poses is that the escape speed from most galaxies is < 1000 km s _1 
(see Figure 2 of iMerritt et al.ll2004l ). and the escape speed from the central bulge is even 
smaller. Therefore, if large recoil speeds are typical, one might expect that many galaxies that 
have undergone major mergers would be without a black hole. This is in clear contradiction 
to the obs ervation that galaxies with bulges all appear to have central supermassive black 



holes (see iFerrarese fe Fordl 120051 ) . It therefore seems that there is astrophysical avoidance 



of the types of supermassive black hole coalescences that would lead to kicks beyond galactic 
escape speeds. From the numerical relativity results, this could happen if (1) the spins are 
all small, (2) the mass ratios of coalescing black holes are all much less than unity, or (3) the 
spins tend to align with each other and the orbital angular momentum. 

The low-spin solution is not favored observationally. X-ray observations of several active 
galactic nuclei reveal relatiyistically broadened Fe Ka fluorescence lines indicative of spins 
a > .9 ( Iwasawa et al.|[l996 : Fabian et al. 2002 ; Reynolds fe Nowak 2003 ; Brenneman fe Reynolds 
20061 ). A similarly broad line is seen in the stacked spect ra of active galactic nuclei in a long 
exposure of the Lockman Hole (jStreblyanska et al.ll2005l ). More generally, the inferred aver- 
age radiation efficiency of superma ssive black holes suggests that they tend to rotate rapidly 



( ISoltanlll982l ; lYu fe Tremaindl2002l ; see lMarconi et al.ll2004l for a discussion o f uncertainties) 



This is also consistent with predictions from hierarchical merger models (e.g. 
2005|). 



Volonteri et al. 



Mass ratios much less than unity may occur in some me rgers, and if the m asses are 
dif ferent enough then the kic k speed can be small. For example, IBaker et al.l (120071 ). followed 
by ICampanelli et al.l (j2007bl ). suggest that the spin kick component scales with mass ratio 
q = mi/m 2 < 1 as q 2 (a2 — qa\)/(l + q) 5 , hence for a\ = — a 2 = 1 the maximum kick speed 
is oc q 2 /(l + q) A . For q < 0.1 this scales roughly as q 2 and hence kicks are small. However, 
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for q > 0.2 the maximum kick is within a factor ~ 3 of the kick possible for q — 1. An 
unlikely conspiracy would thus seem to be required for the masses always to be different by 
the required factor of several. Some tens of percent of galaxies appear to have undergone at 
least one merger with mass r atio > 0.25 within redshift z < 1 (fo r recent observational results 
wit h different method s, see iBell et al.l l2006bl ; lLotz et al.l 120061 . and for a recent simulation 



see iMaller et al.ll2006l ). The well-established tight correlations b etween central black h ole 
mass and galactic properties such as bulge velocity dispersion (see iFerrarese fc Fordll2005l for 
a review) then suggest strongly that coalescence of comparable-mass black holes should be 
common. 

The most likely solution therefore seems to be that astrophysical processes tend to align 
the spins of supermassive black holes with the orbital axis. This astrophysical alignment 
is the subject of this Letter. Here we show that gas-rich mergers tend to lead to strong 
alignment of the spin axes with the orbital angular momentum and thus to kick speeds 
much less than the escape speeds of sizeable galaxies. In contrast, gas-poor mergers show 
no net tendency for alignment, assuming an initially uniform distribution of spin and orbital 
angular momentum vectors. We demonstrate this aspect of gas-poor mergers in § 2. In 
§ 3 we discuss gas-rich mergers, and show that observations and simulations of nuclear 
gas in galactic mergers suggest that the black holes will be aligned efficiently. We discuss 
consequences and predictions of this alignment in § 4. 



2. Gas-poor Mergers 



Several recent models and observations have been proposed as evidence that some ffalac- 
tic mergers occur without a significant influence of gas . Pos sible signatures include the 



Naab k, Ostrikerl 120071 ) and slow rotation and the 



metal richness of giant ellipticals (e.g., 
presence of boxy orbits in th e centers of some elliptical galaxies (e.g., 
Naab. Jesseit. fc Burkertll2006h . 



Bell et al.l l2006al : 



Consider such a gas-free merger, and assume that we can therefore treat the gradual 
i nspir al of two spinning black holes as an isolated system. As laid out clearly by ISchnittman 
(I2004J ). throughout almost the entire inspiral the r e is a strong hierarchy of time scales, 
such that ti nS pi ra i t pr ccess ^orbit- ISchnittmanl (120041 ) therefore derived orbit-averaged 
equations for the spin evolution in the presence of adiabatic dissipation. Such effects can 
lead to relaxation onto favored orientations. The question is then whether, with the uniform 
distribution of orbital and spin directions that seems expected in galactic mergers, there is 
a tendency to align in such a way that the net kicks are small. 
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Using equations A8 and A10 from ISchnittmanl (120041 ). we have evolved the angles be- 
tween the two spin vectors, and between the spins and the orbital angular momentum. We 
find that for isotropically distributed initial spins and orbits, the spins and orbits at close 
separati on are also c l ose to isotropically distributed (see Figured]). Thus, although (as we 
confirm) ISchnittmanl (120041 ) showed that for special orientations the spins might align (e.g., 
for an initial cos#i ~ 1, or as we also discovered, for an initial cos 82 ~ —1), the initial 
conditions resulting in such alignment are special and subtend only a small solid angle. 

The conclusion is that gas-poor mergers alone canno t align spins sufficiently t o avoid 
large kicks due to gravitational radiation recoil. Indeed, ISchnittman fc Buonannol (120071 ) 
find that for mass ratios q > 0.25, spin magnitudes di = d 2 = 0.9, and isotropic spin 
directions, ~ 8% of coalescences result in kick speeds > 1000 km s" 1 and ~ 30% y ield speeds 
> 500 km B _1 . The hieh maximum speeds inferred by ICampanelli et al.l (l2007bi ) are likely 
to increase these numbers. We now discuss gas-rich mergers, which can naturally reduce the 
kick speeds by aligning black hole spins with their orbital axis. 



3. Gas-rich Mergers 



Consider now a gas rich environment, which is common in many galactic mergers. The 
key new element is that gas accretion can exert torques that change the direction but not the 
magnitude of the spin of a black ho le, and that the lever arm for these torques can be tens 



of thousands of gravitationa l radii dBardeen fc Petterson 



(If999h 



1975 



Nataraian fc Armitage 1999 ). In particular, Natarajan Sz Pringle fjl998 ) and Natarajan fc Armitage 



Nataraian fc Pringle 



1998 



demonstrate that the black hole can align with the larger scale accretion disk on a 
timescale that is as short as 1% of the accretion t ime. An important ingredient of this 
scenario is the realization by iPapaloizou fc Pringld (119831 ) that the warps are transmitted 
through the disk on a timescale that is shorter by a factor of l/2a 2 compared with the trans- 
port of the orbital a ngular momentum in flat d isks, where a ~ 0.01 — 0.1 is the standard 
viscosity parameter (jShakura fc Sunyaevl Il973l ). The question that distinguishes gas-rich 
from gas-poor mergers is therefore whether the accreted mass is ~ 0.01 — O.lMbh during the 
sinking of the black holes towards the center of the merged galaxy, where Mbh is a black hole 
mass. 

Numerical simulations show that galactic mergers trigger large gas inflows into the 
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find that the strong gas inflows observed in cooling and star formation simulations always 
produce a rotationally supported nuclear disk of size ~ 1 — 2 kpc with peak rotational 
velocities in the range of 250—300 kins -1 . 

The results of numerical simulations are in good agreement with observations, which 
also show that the total mass of the gas accumulated in the central region of merger galaxies 
can reac h 10 9 — 10 10 M M a nd in some cases account for about half of the enclosed dynamical 



mass (ITacconi et al.lll999l ). Observations imply that the cold, molecular gas settles into a 
geometrically thick, rotating structure with velocity gradi ents similar to these obta ined in 
simulations and with densities in the range 10 2 — 10 5 cm -3 (IDownes fc Solomonlll998l ). Both 
observations and simulations of multiphase interstellar matter with stellar feedback show a 
broad range of g as temperatures, where the la rgest fraction of gas by mass has a temperature 
of about 100 K JvVada fc NormanlboOliliooi l. 



We therefore consider an idealized model based on these observations and simulations. 
In our model, the two black holes are displaced from the center embedded within the galactic- 
scale gas disk. We are mainly concerned with the phase in which the holes are separated by 
hundreds of parsecs, hence the enclosed gas and stellar mass greatly exceeds the black hole 
masses and we c an assume that the black h oles interact i ndepe ndent ly with the disk. Based 
on the results of lEscala etaD tooA hopj ). iMaver et all J2006h . and botti. Colpi. fc Haardt 
(120061 ) the time for the black holes to sink from these separations to the center of the disk 
due to dynamical friction against gaseous and s tellar backgr ound is < 5 x 10 7 yr, which is 
comparable to the starburst timescale, ~ 10 8 yr (jLarsonl 119871 ). 



The accretion onto the holes is mediated by their nuclear accretion disks fed from the 
galactic scal e gas disk at the Bondi rate, MBondi, as long as it does not exceed the Eddington 
rate, M Edd (IGoodman fc Tanl 12004 ). Locally, one can estimate the Bondi radius i?Bondi 



GM hh /Vg » 40 pc (M bh /10 8 M Q ) (vg/100 km s -1 ) -2 that would be appropriate for a total gas 
speed at infinity, relative to a black hole, of v g (we use a relatively large scaling of 100 km s _1 
for this quantity to be conservative and to include random motions of gas clouds as well as 
the small thermal speed within each cloud). The accretion rate onto the holes will then be 
min(M Bondi ,M Edd ), where M Bondi « 1 M Q yr" 1 (t; g /lOOkms- 1 )~ 3 (n/lOOcm- 3 )(M bh /lO 8 M ) 2 . 
We also note that clearing of a gap requires accretion of enough gas to align the holes with 
the large-scale gas flow. 

At this rate, the holes will acquire 1-10% of their mass in a time short compared to the 
time needed for the holes to spiral in towards the center or the time for a starburst to deplete 
the supply of gas. The gas has significant angular momentum relative to the black holes: 
analogous simulations in a planetary formation enviro nment suggest that the circular ization 



radius is some hundredths of the capture radius (e.g.. lHamilton fc Burnslll99ll . Il992l ). This 
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corresponds to more than 10 5 gravitational radii, hence alignment of the black hole spin axes 
is efficient (and not antialignment, since the cumulative angular momentu m of the accretion 



disk is much greater th an the angular momentum of the black holes; see iKing et al.l 12005 
Lodato fc Pringle!l2006h . 



If the black hole spins have not been aligned by the time their Bondi radii overlap and a 
hol e is produc ed in the disk, further alignment seems unlikely (for a different interpretation 
200J). The reason is twofold: the shrin king of the binary due t o circumbinary 



sec 



Liu 



torques is likely to occur within < few x 10 7 yr (jEscala et al.l I2004J . 120051 ). and accretion 
across the gap only occurs at ~ 10% of the rate it would have for a single black hole 
( Lubow. Seibert. fc Artymowicz 1999 ; Lubow fc D'Angeloll2006 ; MacFadyen &: Milosavljevic 



20061 ). with possibly even smaller rates onto the holes themselves. This therefore leads to the 
gas-poor merger scenario, suggesting that massive ellipticals or ellipticals with slow rotation 
or boxy orbits have a several percent chance of having ejected their merged black holes but 
that other galaxy types will retain their holes securely. 



4. Predictions, Discussion and Conclusions 



We propose that when two black hole accrete at least ~ 1 — 10% of their masses during 
a gas-rich galactic merger, their spins will align with the orbital axis and hence the ultimate 
gravitational radiation recoil will be < 200 km s _1 . In this section we discuss several other 
observational predictions that follow from this scenario. The best diagnostic of black hole 
spin orientation is obtained by examining AGN jets. All viable jet formation mechanisms 
result in a jet that is initially launched along the spin axis of the black hole. This is the 
case even if the jet is energized by the accretion disk rather than the black hole spin since 
the orientation of the inn er accretion disk will be sla ved to the black hole spin axis by the 
Bardeen-Petterson effect (IBardeen fc Pettersonlll975l ). 



At first glance, alignment of black hole spin with the large scale angular momentum of 
the gas would seem to run contrary to the observation that Seyfert galaxie s have jets that are 



randomly oriented relative to their host galaxy disks ( Kinney et al 



morphology is not consistent with recent major mergers (jVeilleux 



2000). However, Seyfert 



20031 ). hence randomly- 



oriented minor mergers or internal processes (e.g., scattering of a giant molecular cloud into 
the black hole loss cone) are likely the cause of the current jet directions in Seyfert galaxies. 

Within our scenario, one will never witness dramatic spin orientation changes during 
the final phase of black hole coalescence following a gas-rich merger. There is a class of 
radio-loud AGN known as "X-shaped radio galaxies", however, that possess morphologies 
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interpreted precisely as a rapid (< 10 5 yr) re-alignment of black hole spin during a binary 



black hole coalescence (lEkers et al. 
20031 : lKomossal[2003bl : lLal fc Raol 



1978 



2007 



Dennett-Thorpe et al.l 120021 ; IWang. Zhou k, Dong 



Cheungil2007l ). These sources have relatively normal 



"active" radio lobes (often displaying jets and hot-spots) but, in addition, have distinct 
"wings" at a different position angle. The spin realignment hypothesis argues that the wings 
are old radio lobes associated with jets from the one of the pre-coalescence black holes in 
which the s pin axis possessed an en tirely different orientation to the post-coalescence remnant 
black hole (IMerritt fc Ekerd 120021 ) . If this hypothesis is confirmed by, for example, catching 
one of these systems in the small window of time in which both sets of radio lobes have 
active hot spots, it would contradict our scenario unless it can be demonstrated that all such 
X-shaped radio-galaxies originate from gas-poor mergers. 

However, the existence of a viable alternative mechanism currently prevents a compelling 
case from being made that X-shaped radio sources are a unique signature of mis-aligned 
black hole coalescences. The collision and subsequent la teral expansion of t he radio galaxy 
backflows can equally well produce the observed wings (ICapetti et al.l 120021) . Ind eed, there 



is circumstantial evidence supporting the backflow hypothesis. iKraft et al.l (120051 ) present a 
Chandra observation of the X-shaped radio galaxy 3C 403 and find that the hot ISM of the 
host galaxy is strongly elliptical, with a (projected) eccentricity of e ~ 0.6. Furthermore, the 
wings of the "X" are closely aligned with the minor axis of the gas distribution, supporting a 
model in which the wings correspond to a colliding backflow that has "blown" out of the ISM 
along the direction of least resistance. Although 3C 403 is the o nly X-shaped radio g alaxy 
for which high-resolution X-ray maps of the hot ISM are available, ICapetti et al.l (120021 ) have 
noted that a number of X-shaped sources have wings that are oriented along the minor axis 
of the optical host galaxy. This suggests that the conclusion of IKraft et al.l (120051 ) for 3C 403 
may be more generally true. 

There is one particular system, 0402+379, that might provide a direct view of spin 
alignment in a binary black hole sys tem. Very Long Baseline Array (VLBA) imaging of this 
radio galaxy by iManess et al.l (120041 ) discov ered two compact flat sp ectrum radio cores, and 
follow-up VLBA observations presented by iRodriguez et al.l (120061 ) showed the cores to be 
stationary. A binary supermassive black hole is the most satisfactory explanation for this 
source, with the projected distance between the two black holes being only 7.3 pc. Within 
the context of our gas-rich merger scenario, these two black holes already have aligned spins. 
Existing VLBA data only show a jet associated with one of the radio cores. We predict 
that, if a jet is eventually found associated with the other radio core, it will have the same 
position angle as the existing jet. 



Evidence already exists for alignment of a coalescence remnant with its galactic scale 
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gas disk. iPerlman et al.l (120011 ) imaged the host galaxies of three compact symmetric objects 
and discovered nuclear gas disks approximately normal to the jet axis. The presence of such 
a nuclear gas disk as well as disturbances in the outer isophotes of all three host galaxies 
suggests that these galaxies had indeed suffered major gas-rich mergers within the past 
10 8 yr. 

In conclusion, we propose that in the majority of galactic mergers, torques from gas 
accretion align the spins of supermassive black holes and their orbital axis with large-scale 
gas disks. This scenario helps explain the ubiquity of black holes in galaxies despite the 
potentially large kicks from gravitational radiation recoil. Further observations, particularly 
of galaxy mergers that do not involve significant amounts of gas, will test our predictions 
and may point to a class of large galaxies without central black holes. 
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UMCP-Astronomy Center for Theory and Computation Prize Fellowship program for sup- 
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Fig. 1. — Demonstration that in gas-poor mergers there is no tendency to align black hole 
spin angles. Here we show the distribution of the dot products with the orbital angular 
momentum axis of the final spin axis of the larger mass (cosfli and smaller mass (cos#2,/i) 
black hole, evolved using the formalism of ISchnittmanl (120041 ) . We assume that at an initial 
separation of 1000m (where m = 1 is the total mass of the binary) the spin directions 
and orbital axis are distributed isotropically, and we integrate inward to 10m assuming 
component masses m\ = 0.55 and m 2 = 0.45 and dimensionless spin parameters di = a,2 = 1. 
The final spin angles show no alignment towards each other or towards the orbital axis, hence 
other mechanisms are needed to avoid ejecting merged black holes from their host galaxies. 



